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We discuss symmetry breaking quantum phase transitions on the oft studied Bethe lattice in 
the context of the ferromagnetic scalar spherical model or, equivalent ly, the infinite Nf limit of 
ferromagnetic models with 0{Nf) symmetry. We show that the approach to quantum criticality is 
characterized by the vanishing of a gap to just the global modes so that all local correlation functions 
continue to exhibit massive behavior. This behavior persists into the broken symmetry phase even 
as the order parameter develops an expectation value and thus there are no massless Goldstone 
bosons in the spectrum. We relate this feature to a spectral property of the graph Laplacian shared 
by the set of "expander" graphs, and argue that our results apply to symmetry breaking transitions 
on such graphs quite generally. 



I. INTRODUCTION 

The study of statistical mechanical systems on the 
Bethe lattice dates back to Bethe's early work on binary 
alloys in the 1930&^. Formally, the Bethe lattice is a reg- 
ular graph with no loops in the infinite size limit; its tree- 
like local structure makes it amenable to real space self- 
consistent treatments for systems with short ranged in- 
teractions without entirely losing the notion of distance, 
as happens if instead we resort to infinite ranged interac- 
tions. Specifically, for ferromagnetic, and hence unfrus- 
trated, models it makes the Bethe-Peierls approximation 
exact. For frustrated systems, such as spin glasses, the 
simplification achieved is more problematic and has been 
the subject of much ongoing work^. Altogether, there is 
a large set of results on various aspects of classical sta- 
tistical mechanics on the Bethe lattice. 

Quantum statistical mechanics on the Bethe lattice has 
received less attention by comparison. One early set of 
results concern the Hubbard model on the Bethe lattice 
in the infinite coordination limit^-. More recently, purely 
bosonic models on the Bethe lattice have been the subject 
of investigatio n wi th the introduction of the quantum 
cavity met hod. f^*^ This has led to computational results 
on the quantum Ising model, the quantum Ising spin glass 
and on the Bose-Hubbard model. In addition, we note 
variational work on the ground state of the quantum Ising 
modeF and on a spherical model of a spin glass^. 

In this paper we study the quantum unfrustrated (fer- 
romagnetic) problem on the Bethe lattice analytically. 
For technical reasons, we do this nominally in the form of 
the spherical model for a scalar field but it is essentially 
also the large Nf limit of the nearest neighbor 0{Nf) 
quantum rotor model on the Bethe lattice. We find some 
surprising results from the perspective of quantum phase 
transitions on Euclidean lattices. First, we find that the 
critical point is marked by a single global mode (or Nf 
modes in the 0{Nf) interpretation) descending in en- 
ergy to become degenerate with the ground state while 
all other states remain a finite distance away in energy. 
Second, no local operators couple to this mode in the 
thermodynamic limit so that all local correlation func- 
tions exhibit a gap at criticality. Alternatively, once the 



thermodynamic limit has been taken, the response to a 
field on any finite subvolume is bounded even if that 
volume is taken to infinity. Third, this state of affairs 
persists into the broken symmetry phase in that while 
there are 2 (or 0{Nf) worth of) degenerate broken sym- 
metry states as required by the global symmetry, there 
are no Goldstone bosons and all local correlators remain 
gapped. We believe our results are robust to moving 
away from the spherical/ A^/ = oo limit. Indeed, the 
Nf = 2 problem is the particle-hole symmetric transi- 
tion in the Bose-Hubbard model studied using the cavity 
method by Semerjian et aP whose results imply that a 
macroscopic superfluid density develops at the quantum 
phase transition (QPT) without closing the excitation 
gap in the single site correlation function. 

These results should surprise the reader. On the face of 
it, the lack of Goldstone behavior confounds the standard 
intuition regarding broken symmetry phases. Aficiona- 
dos of classical models on the Bethe lattice may counter, 
rightly, that their phase transitions are accompanied by 
exponential rather than power law correlations, as this is 
sufficient for the classical susceptibility to diverge. How- 
ever, if one computes the quantum susceptibility - by 
which we mean the response to an infinitesimal field ap- 
plied over a subvolume whose size is taken to infinity after 
the thermodynamic limit - it remains finite, because the 
matrix elements in question do not couple to the global 
mode. The analogous susceptibility of the classical model 
diverges to signal the phase transition. 

We turn now to a brief discussion of what we mean 
by the Bethe lattice as it plays an important role in the 
following. The traditional, and simplest, way to define 
the Bethe lattice is through a sequence of Cayley trees. 
The Bethe lattice then corresponds to the interior of a 
large tree where each site is connected to exactly z neigh- 
bors. This construction makes clear that the number of 
sites within a fixed distance of a given site scales expo- 
nentially with distance and thus connects, qualitatively, 
with the infinite dimensional limit of hypercubic lattices 
that is also often invoked in the context of making self- 
consistent treatments exact. Unfortunately, the bound- 
ary of a finite Cayley tree is always a finite fraction of 
the bulk whence the choice of boundary conditions can 
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complicate the thermodynamic hmilF. Indeed, the differ- 
ence between a ferromagnet and a glass can be phrased 
solely in terms of boundary conditions. To get around 
this problem, one can proceed differently and consider 
instead the ensemble of regular graphs with N nodes. 
This consists of all graphs where each site is connected 
to exactly z nearest neighbors. Picking a graph at ran- 
dom from this ensemble in the N ^ oo limit yields an 
alternative definition of the Bethe lattice which is known 
to be perfectly satisfactory for the classical ferromagnetic 
problem^SI that is what we shall use in this paper. 
At finite members of the ensemble contain loops of 
characteristic size logA^ and the graphs are not entirely 
homogeneous. However, homogeneity is restored in the 
limit ^ oo as the loop size diverges. 

With this definition in hand, we can comment on the 
feature of the Bethe lattice that brings about the surpris- 
ing features we report: the boundary of any subvolume of 
the lattice is always a finite fraction of the whole whence 
(arbitrarily) localized excitations cannot be made low en- 
ergy. In more formal terms, the Bethe lattice has what 
is known to graph theorists as a positive Cheeger con- 
stant h - that is the minimal ratio of boundary to bulk 
of its nontrivial subgraphs. A theorem - the Cheeger 
bound - guarantees that the Laplacian has a gap of at 
least h'^ /2 between the uniform Perron- Frobenius ground 
state and the rest of the spectrum. Moreover, this prop- 
erty is shared by a large collection of so-called expander 
graphs^^, of much interest in quantum and classical in- 
formation theory, and we expect our qualitative results 
to be shared by models on any such graph. 

In the balance of this paper, we first introduce and 
solve the classical spherical model on the Bethe lattice 
which introduces most of the ingredients needed in a sim- 
pler setting. This classical model on the Bethe lattice has 
been solved previousl}ff^ by a different technique. We 
next introduce and solve the quantum spherical model 
which exhibits the features we described above. (We note 
that the quantum spherical model on Euclidean lattices 
has been solved by Vojta^^.) We end with a discussion 
and comments on the generality of our results. 



II. THE CLASSICAL SPHERICAL MODEL 

The classical spherical model on a z-regular graph Q 
with N nodes is defined by the Hamiltonian 
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(1) 



subject to the global constraint that (j)^ = N. Here, 
the graph Laplacian is given by 



(2) 



A is the adjacency matrix of ^, and the overall shift of z 
ensures that H is positive semidefinite. Alternatively, we 
can view this choice of the Hamiltonian as the one which 



gives a system of coupled oscillators with nonnegative 
frequencies. 



A. Phase transition 

The partition function of the model (at inverse tem- 
perature /3 = 1/T) is given by the expression 



(3) 



Representing the delta function by a Lagrange multiplier 
A and performing the Gaussian integral over the 0^, we 
obtain the effective action 



Z{/3) = J ^;^g-A^[27vTrlog{/3/:-2iAl}+iA] 



(4) 



As A/" ^ oo, this integral may be performed by steepest 
descent. The effective propagator for the field is given 
by (defining -/i/T = 2iX): 

=T(/: + /il)r.i (5) 

and the constraint that = N gives the self- 

consistency condition 

Ar = Tr[T(£ + /il)-i] (6) 

We may rewrite this in terms of the eigenvalues of 
the graph Laplacian 
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(7) 



where we have explicitly written /i(T) to emphasize 
that self-consistency forces /i to depend on temperature. 
Throughout the paper, we use indices a^P to label the 
modes of the Laplacian jC with energy and eigenvector 
uf and i, j to refer to sites of the lattice. 

The spectrum of z-regular random graph Laplacians 
and the related problem of hopping on the Bethe lat- 
tice have been extensively studied^^^^. The spectrum 
consists of a continuum of states in the interval [z — 
2\/z — 1, z + 2\/z — 1], with the density of states 



Pc(e) = 



z v/4(z - 1) - (6 - zy 
2tt z^-{e- zY 



(8) 



However, the above density of states does not include the 
uniform Perron- Frobenius eigenvector, which has ampli- 
tude 1/\/N on each site and energy 0. Thus, the full 
density of states is 



(5(e) + p,{e) 



(9) 



The existence of this spectral gap, and the uniqueness of 
the low-lying state is guaranteed by the fact that the z- 
regular random graphs have a positive Cheeger constant. 
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We now return to the self-consistency equation Eq. |7|) 
and rewrite it in terms of the density of states Eq. ([9f: 



1 
T 



1 



I 



de 



(10) 



-2^(^-1) 



The usual argument for Bose-Einstein condensation now 
follows: as T ^ 0, we must decrease /i in order to satisfy 
this sum rule. At high temperatures, the thermodynamic 
limit can be taken straightforwardly without paying es- 
percial attention to the uniform state, and one can satisfy 
the self-consistency equation. However, since /i cannot 
decrease below the lowest eigenvalue of the Laplacian - 
as this would render the steepest descent calculation un- 
stable - the smallest that fi can be is zero. This occurs 
at a critical temperature Tc given by 



de 



\z-2.yiz-l) 



J 



z{z - 2) 
z-l 



(11) 



This result for the transition temperature was previ- 
ously derived in Ref. 12, Clearly, Tc = for z = 2, and 
there is no finite-temperature transition on the chain, as 
expected. However, Tc is finite for 2: > 2, and below this 
temperature, we cannot satisfy the self-consistent equa- 
tion in the naive thermodynamic limit discussed above. 
Rather, we must keep track of the uniform mode in (10) 



and take 11 = ^ {l - + 0{N-'^) for T < Tc. 

For comparison, in the 0{Nf) model, the Lagrange 
multiplier A would become a field A^, which at finite N 
on the random graph need not be homogeneous at the 
large Nf saddle point. However, the infinite Bet he lattice 
is homogeneous and so in the thermodynamic limit we 
expect to recover the same saddle point as in the simpler 
spherical model. 



B. Ordered phase 

How does the system behave in the low temperature 
phase? Consider the correlation function computed in 
the eigenbasis of £, namely 



TSaf: 



(12) 



Below the critical point, the amplitude of the lowest 
eigenmode is given by 



N 1 



T 



(13) 



The thermal occupation of the uniform state is macro- 
scopic, indicative of long range order in the system. In 
light of this, we revisit the steepest descent calculation. 



but this time treat the lowest mode separately. First, 
integrate out the TV — 1 higher modes 

Z{(3) = J dA#oe'^(^°"^) 



/n< 
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-N 



— / d<j)Qd\e 



ix[ 1 



+ 27V E log{/3e«-2iA} 



N 1^ 2N 



and now perform the A integral by steepest descent, giv- 
ing the self-consistency equation (again ijl/T = —2iX) 



e + /X 



(14) 



For T near Tc, ji <^ z — 2\J z — 1 so we expand in powers 
of /i and solve: 



1 



hT \N 



(15) 



where h = J^^^^ ^ ' 

Using this solution and approximating the steepest de- 
scent value of the integrand to quadratic order in /i((/)o), 
we find that 

Z{(3) ^ J d^oe-^^^ff^^'^ (16) 

where the effective potential for ^0 is (dropping con- 
stants): 



1 



26T2 



^-1 



(17) 

Thus, we see that the effective potential has Landau- 
Ginzburg form, and predicts a mean field symmetry 
breaking for T < Tc such that (^0) ~ V^. That is. 



C. Spatial correlations 



We now consider the spatial behavior of the Green's 
•j). As this is ultimately a Gaussian 



function Gij 



theory, the only interesting correlator is the two-point 
function. Since we do not know the exact eigenvectors 
of the graph Laplacian at finite we cannot obtain Gij 
from Gaf3 directly. Instead, we work directly in the ther- 
modynamic limit of the Bethe lattice and assume that 
Gij = G{\i — j\) is translation invariant. In the disor- 
dered phase, Gij satisfies the equation of motion 



{Cii^ + li8ii')Gi' 



(18) 
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Using the ansatz Gij = Ae ■'1/^, one finds 



2{z-l) 



[z - 2){z + + z.J{z + ijlY -4.{z-l) 

^ {z + ^^)-^{z + ^i)^-A{z-l) \'~'\ 
\ 2{z-l) ) ^^'^ 

which, as /i ^ 0+ at the phase transition, reduces to 
z-1 



z{z-2) 



(z-1) 



(20) 



This agrees with the weh-known correlation length ^ = 
l/ln(z — 1) of the Bethe lattice at criticality and corre- 
sponds to a diverging global susceptibility. 

Since the uniform mode has vanishing weight in the 
disordered phase, it makes no contribution to the Green's 
function for T > Tc. In the low temperature phase, the 
uniform mode has macroscopic occupation correspond- 
ing to a non-zero value of (^^). In this case, the ex- 
pression (12), calculated at /i = 0, corresponds to the 



connected correlation function G^^- = {(t)i(})j) — 
throughout the low temperature phase. That the corre- 
lation length sticks at its critical value and the suscepti- 
bility remains divergent throughout the condensed phase 
is the Bethe lattice analogue of the classical Goldstone 
theorem, which arises quite naturally when this model is 
viewed from the large TV/ point of view. 



III. THE QUANTUM SPHERICAL MODEL 

The quantum spherical model is given by the 
Hamiltonian^^ 



(21) 



with canonical commutation relations [^i,7rj] = i6ij, and 
the mean spherical constraint: 



N 



(22) 



As discussed in Ref. 13, the partition function can be 
rewritten in functional integral form as 

ij \ i / 

where r is an imaginary time parameter. Transforming 
to frequency space and expanding in eigenvectors of £, 



^(^,/5) 



X e 



d\Wd(\)^(uj) (23) 

OL 

\ Y: (/>«(a;)[/3(^+e«)-2iA]0«(-a;)-iAAr 



where we sum over discrete frequencies = 
The self-consistency equation is now given by 



/3 • 



TV 



1 



Performing the frequency summation, we obtain 



(24) 



^4e 



V9 



+ m) (25) 



This result was computed directly within the Hamilto- 
nian formalism in Ref.lT3l Note that as ^ ^ 0, we recover 
the classical result Eq. 



A. Quantum phase transition 



At T = 0, the self-consistency equation reduces to 



1 



N ^ 



1 



de 



(26) 



Following a similar argument as in the classical case, we 
find that the system has a quantum critical point at the 
critical coupling defined by 



-1/2 



/ 



de 



2v^ 



(27) 



This integral is convergent, and can be done numerically 
to find the precise value of g^ Again, for g < gc^^ must 
be careful regarding the thermodynamic limit and take 

V~9 



2N 



(i-v^) ' 



(28) 



in order to determine the macroscopic occupation of the 
zero mode. As expected, the quantum spherical model 
undergoes a Bose condensation transition at the critical 
coupling gc. 



B. Correlations 

The imaginary time Green's function at zero tempera- 
ture is most easily computed by taking a limit from finite 
temperature. From the partition function Eq. (23), we 
find 



f duo e'^^ 



(29) 
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tion Gij{uj) is given by the mode sum 



FIG. 1: Support of Gij{ijo) at zero temperature for z = 3. The 
isolated (blue) line indicates the contribution of the uniform 
mode which vanishes as l/A/" above Qc (dashed) but is finite 
below Qc signaling the long-range order (solid) . Obtained from 
numerical inversion of Eq. (26). 



Each mode a decays exponentially in imaginary time at 
a rate corresponding to the gap to exciting that mode. 
Thus, the spectral response of the time-ordered propaga- 
tor in re al frequen cy space is precisely a delta function 
at = v^^(e + /i). 

Now we are in a position to understand the peculiar- 
ity of the quantum condensation transition on the Bethe 
lattice. The single site Green's function is given by trans- 
forming Gcx[3 back to position space using the eigenmodes 
of the graph Laplacian: 



1 ^ ^( 



V9 



(30) 



a>0 



- \/g{ea ^ IJ)) 



Without detailed knowledge of the eigenvectors of the 
Laplacian but merely its density of states, this formula 
already allows us to sketch the spectral response of local 
single particle excitations as in Fig. [l] As expected, the 
quantum phase transition is signaled by the continuous 
closing of a spectral gap to a uniform mode. The unusual 
feature is that this state is isolated from the remainder 
of the spectrum and its weight vanishes as 1/N in the 
disordered phase. Meanwhile, in the ordered phase, the 
weight 



1 Vfl 

N2^ 



(31) 



is finite, reflecting the condensation into the uniform 
state. 

The general spectral features of correlations sketched 
above may be supplemented in detail by assuming ho- 
mogeneity of the infinite Bethe lattice and exploiting 
our knowledge of the classical spherical model. The fre- 
quency resolution of the zero temperature Green's func- 



Gij{uj) 



ufuf 



Oi g 



(32) 



which has the same form as the mode sum for a classi- 

2 

cal model at chemical potential /i + — . Thus, we use 



the classical Green's function in Eq. (19) and make this 
substitution for ji. The formal expression is rather un- 
enlightening but there is much information in the pole 
structure (see Fig. [2|. The most important feature is a 
pair of vertical square root branch cuts on the imaginary 
axis corresponding to the bulk spectrum. These never 
pinch off the real axis as g is varied; rather, the closest 
branch point sets the dominant decay rate, which agrees 
exactly with the previous discussion. 

At finite N the branch cuts break into lines of poles and 
in addition one should recall the contribution of the uni- 
form mode. This provides a pair of isolated simple poles 
at ±i^/gjl with residue in the disordered phase. As 
g ^ gt 1 these poles pinch off the real axis corresponding 
to the closing of the gap and the phase transition. In the 
condensed phase, these poles have non- vanishing residue 
indicating the long range order in imaginary time. 

Finally, we note that the static susceptibility x of the 
quantum model to a global field (applied before the ther- 
modynamic limit) is simply given by the zero frequency 
lattice sum of Gij{uo)^ which is equivalent to summing 
Eq. (19) over the lattice. This susceptibility diverges 
throughout the broken symmetry phase, as usual. How- 
ever, taking a large subvolume and computing its suscep- 
tibility after the thermodynamic limit gives a convergent 
result, reflecting the absence of Goldstone bosons in the 
spectrum. 



IV. DISCUSSION 

There are two salient qualitative features of the spher- 
ical model condensation transition on the Bethe lattice 
that we have derived above: (1) despite its vanishing 
weight in the response to local excitations, there is a 
global mode closing the many-body gap in the disordered 
phase. Alternatively, this can be viewed as a 'midgap' 
mode driving the transition, where by midgap we mean 
in the middle of the gap to local excitations. (2) It is 
via macroscopic occupation of this mode that long range 
order develops in the symmetry-broken phase. However, 
local excitations remain gapped. In fact, not only are 
there no Goldstone bosons, the spectral response of any 
operator with bounded support will remain gapped. 

We believe these features are quite general to unfrus- 
trated transitions on the Bethe lattice and expander 
graphs. For one thing, as already noted in the intro- 
duction, our results are entirely consistent with those of 
Ref. 5 on the Bose-Hubbard model wherein the particle- 
hole symmetric transition is the 0{Nf) model at Nf = 2. 
This indicates that the features we have found at Nf = oo 
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FIG. 2: Analytic structure of Gij{<jj) in the complex plane. 
Wiggly lines (green) indicate branch cuts. The two crosses 
(red) indicate the poles due to the zero mode with residue 
1/N. The following formulae for cut locations are given 
at /i = (criticality): The hyperbolic cuts are given by 
Re[a;]^ — Im[a;]^ = —zg. The vertical cuts extend between 

g{z db 2^z — T). These cuts control the slowest decay 
since we can deform the contour past the uniform mode pole 
in the thermodynamic limit. 

are robust throughout the range of 0{Nf) models and 
not merely artifacts of the spherical model limit. We be- 
lieve this claim will be susceptible to a proof in the 1/Nf 
expansion^^. Furthermore, the results should apply to 
sufficiently weakly disordered ferromagnetic models, for 
which the disorder does not close the mode gap. 

More intuitively, these results should be stable to the 
inclusion of self-interactions of the scalar field beyond the 
self-consistent Gaussian theory treated here. First, at the 
critical point, weak interactions cannot close the gap be- 
tween the soft global modes and other massive modes; in- 
deed, they typically cause the low energy density of states 
to decrease even further. Second, the usual heuristic ar- 
gument for a Goldstone mode is seen to lead to a gap 
as follows: consider the simplest case of an 0(2) = U{1) 
broken symmetry in the broken phase. In the usual fash- 
ion we derive an effective Lagrangian for the phase fluc- 



tuations in the symmetry-broken phase by freezing out 
modulus fluctuations in |(^)|: 

^eff « ^ I (0) P E + ^ K^) I' E ^^^V^^ (33) 

This is a quadratic action for a massless free field ^, which 
would lead to gapless behavior of local excitations were 
the Laplacian gapless. However, the Laplacian spectrum 
on expander graphs is gapped to all modulated excita- 
tions. While the global rotation is still, necessarily, of 
zero energy, it no longer follows that there exist local ex- 
citations of arbitrarily low energy. This corresponds to a 
breakdown of our intuitive understanding of generalized 
rigidity. Quite generally then, we conjecture that there 
is an ' ant i- Goldstone' theorem that applies to symmetry 
breaking transitions on expander graphs. 

We should also note that the fully connected models 
for which mean field theory is exact exhibit similarly 
gapped behavior, even after the usual rescaling for ex- 
tensivity. Consider the fully connected Heisenberg ferro- 
magnet whose Hamiltonian is H = —J/N^.j Si ■ Sj = 
— J/NS^^^. Trivially, the gap to the first excited state is 
J and does not vanish in the thermodynamic limit. The 
fully connected graph does not have bounded coordina- 
tion so is not technically an expander graph, however the 
surface to volume ratio is, nonetheless, very large. 

Clearly there are many subtle issues that may invali- 
date the usual physics lore when models are considered 
on non-Euclidean graphs. As quantum models are stud- 
ied more extensively in the context of quantum complex- 
ity theory and on random networks, an enhanced under- 
standing of universal symmetry breaking properties will 
potentially play an important role in their analysis. 
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